Polymer-mediated entropic forces between scale-free objects 
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The number of configurations of a polymer is reduced in the vicinity of a barrier, and the resulting 
loss of entropy leads to a repulsive force. When the obstructions are scale invariant shapes (such 
as cones, wedges, lines or planes) the only relevant length scales are the polymer size R$ and 
characteristic separations, severely constraining the functional form of entropic forces. Specifically, 
we consider a polymer (single strand or star) attached to the tip of a cone, at a separation h from a 
surface (or another cone). At close proximity, such that h <C Ro, separation is the only remaining 
relevant scale and the entropic force must take the form F = AksT/h. The amplitude A is universal, 
and can be related to exponents n governing the anomalous scaling of polymer correlations in the 
presence of obstacles. We use analytical, numerical and e-expansion techniques to compute the 
exponent r\ for a polymer attached to the tip of the cone (with or without an additional plate or 
cone) for ideal and self-avoiding polymers. The entropic force is of the order of O.lpN at 0.1/mi for 
a single polymer, and can be increased for a star polymer. 

PACS numbers: 64.60.F- 82.35.Lr 05.40.Fb 
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I. INTRODUCTION 

A host of single molecule manipulation techniques 0- 
employing atomic force microscopes (AFMs) [(|, mi- 
croneedles 0, optical [H Q and magnetic [l(| tweezers, 
have provided quite detailed studies of shapes and forces 
in lon g po lymers. The positional accuracy of an AFM 
tip (ol lll| can be as good as a few nm, while forces of 
order of 1 pN can be measured, with measurements car- 
ried out in nearly biological conditions [HI, [l3j]. These 
enhanced sensitivities approach the range where entropic 
forces of long polymers in a solvent can be significant even 
when the deformation of the polymer is relatively slight. 
While the main thrust of the experimental research is ex- 
traction of specific information about molecular confor- 
mations from force-displacement behaviors of the poly- 
mers, it is important to take note of polymeric features 
that are independent of the microscopic details [3] . 

At a finite temperature T, a long homo-polymer in a 
good solvent fluctuates between a large number of con- 
figurations, and its 'macroscopic' properties (such as size 
and elasticity beyond a characteristic persistence length) 
are governed by entropic considerations. To illustrate 
entropic elasticity, consider the hypothetical situation of 
a polymer with one end fixed to the origin, while the 
other end explores allowed positions at r. In equilib- 
rium, due to rotational symmetry (r ) — 0. The symme- 
try is broken if a weak force F is applied by a point-like 
probe at the end of the polymer. Then, in ci-dimensional 
space (f ) = {Rl/dk B T)F, where R^ is the mean squared 
end-to-end distance in the absence of external force [3] . 
This is mostly a realignment, without substantial change 
of polymer conformations for F <C F CI = clksT/R^. 
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For i? =0.1/im at room temperature F CT ~ O.lpN is the 
crossover value that separates weak and strong forces; 
stronger force (F ^> F CT ) leads to substantial stretching 
of the polymer and associated reduction in entropy. Sim- 
ilar considerations come into play when the attachment 
points are replaced by macroscopic obstacles (probes): 
Under a strong stretching force polymer configurations 
are not very much influenced by the shapes of the probes, 
and convenient models (such as Ref. }15j) have been de- 
veloped for the resulting force-extension relations. Here, 
we focus on the opposite limit of strong compression in 
which the loss of entropy, and the corresponding force, 
are dominated by confinement within boundaries. 

In this work we consider setups where a polymer is at- 
tached by one end to an object (the probe), while the 
other end is free. The probe is then brought to the vicin- 
ity of a second solid obstacle, and the loss of polymer 
entropy leads to a force between the two objects. In 
particular, we consider confining (probe and obstacle) 
shapes that are self-similar (at least on scales compara- 
ble to Ro) at separations /iCi? . Possible setups include 
sharp AFM tips approaching flat surfaces covered with 
polymers, as in experimental studies (l6l - [20l | of polymer 
brushes. Entropic considerations in a dense brush must 
account for interactions between polymers t 2l|, com- 
plicating the calculation of forces exerted on the AFM 
tip [iH ■ Here we consider the simpler cases of a single lin- 
ear or star polymer, avoiding the dense limit. For h -C Ro 
dimensional considerations suggest the particularly sim- 
ple force-separation relation F = AksT/h, where the 
'universal' amplitude A only depends on basic geometri- 
cal properties, and gross features of the polymer. 

In a previous paper we outlined some results pertain- 
ing to entropic interactions, and critical exponents, for 
a polymer attached to a cone (23[. In this manuscript 
we expand on the details of the calculations, as well as 
generalize and expand the results to other setups. Sec- 



tion |TT] demonstrates how the force amplitude A is re- 
lated to the exponent rj, characterizing the anomalous 
scaling of polymer correlations. The latter exponent de- 
pends on the shape of confining boundaries, and is the 
focus of computations that follow. In Sec. IIIII we outline 
the general approach to computing r\ for ideal 'polymers 
near scale-invariant surfaces, while in Sec. IIVI we make 
explicit computations for cones and plates. The case of 
self-avoiding linear polymers in d = 3 is treated by Monte 
Carlo simulations in Sec. [V] Perturbative e = 4 — d ex- 
pansions are employed in Sec. IVII to compute 77 for linear 
and star self-avoiding polymers attached to cone tips, or 
to a contact point between a cone and a plate. 




II. POLYMER-MEDIATED FORCES 

Let us first consider an idealized setup in which a poly- 
mer is attached to the tip of a solid cone, approaching a 
solid plate (or another cone) . This exemplifies a geometry 
of obstacles in which the only (non-microscopic) length 
scale is provided by the tip-plate (or tip-tip) separation 
h. The polymer itself undergoes self-similar fluctuations, 
spanning length scales intermediate between microscopic 
(monomer size, or persistence length a) and macroscopic. 
The latter is set by the typical end-to-end distance Rq, or 
by the radius of gyration R g (which differs from Rq only 
by a multiplicative factor of order of unity). The typical 
size of the polymer grows with the number of monomers 
TV through the scaling relation Rq ~ R g oc N". The 
self-similar shape of the obstacles also presumably ex- 
tends only up to a characteristic scale H, say the height 
beyond which the cone is terminated or changes shape. 

Neglecting any structure (and hence energy) associated 
with the polymer, variations in free energy as the cone 
(plus tip-attached polymer) approaches the plate are en- 
tirely entropic in origin and proportional to fc^T. The 
entropic force has dimensions of energy divided by length, 
and at separations a <C h <C Rq <C H the relevant length 
scale is h, and it is reasonable to posit 

F = A k -f. (1) 

This is because the polymer configurations are con- 
strained only on the scale of confinement h; increasing 
the length of the polymer or the size of the cone (as long 
as h <C Ro -C H) should not influence the force; further- 
more the entropy change is independent of a. 

Indeed, the simple force law of Eq. (|T|) applies to all 
circumstances where the separation provides the only rel- 
evant length scale, and follows easily from various poly- 
mer scaling forms (see, e.g. the derivation below) such 
as in Refs. [24l426|. The dimensionless amplitude A will 
depend on geometric factors characterizing the confining 
boundaries such as the opening angle of a cone 6 (and if 
tilted, on the corresponding angle). It will also depend 
on factors that characterize the scaling of polymeric fluc- 
tuations, thus differing in cases of ideal and self- avoiding 



(d) (e) 

FIG. 1. (a) Polymer attached to the tip of a solid cone with 
apex semi-angle Q (configuration "c"); positions are described 
by the spherical coordinates r, 8 and azimuthal angle <j> (not 
shown), (b) The tip, where the polymer is attached, is at 
a distance h <C -Ro from the plate, (c) The tip touching 
the plate (configuration "cp"). (d) Tip is at a finite distance 
from a plate to which the polymer is attached, (e) Polymer 
attached to both surfaces. 



polymers, and for linear, star, and brush polymers. In 
the following we shall demonstrate that in a number of 
setups, the amplitude A can be related to variations of 
the (universal, but shape dependent) exponent 77, char- 
acterizing polymeric correlations. 

A chain starting at the tip of a cone provides a first 
approximation to a polymer linked to the tip of an AFM 
probe as depicted in Fig. [TJi. With the cone far away 
from a plate (h 3> Ro), the number of configurations of 
the polymer grows with the number of monomers as 

K = bz N N^^-\ (2) 

where the effective coordination number z, as well as the 
pre-factor 6, depend on the microscopic details (such as 
the scale a), while the 'universal' exponent 7 C only de- 
pends on the cone angle. When the cone touches the 
plate as in Fig. QJ, the number of configurations is re- 
duced to Af cp with the same form as Eq. but with 
a different exponent 7 cp (0). We shall henceforth use an 
exponent subscript 's' (as in 7 S ) to refer to the above 
cases, with "s=c" for cone and "s=cp" for cone+plate; 
the absence of a subscript (as in 7) will signify a free 
polymer. The work done against the entropic force in 
bringing in the tip from afar to contact the plate can 
now be computed from Eq. ([1]) as 

W= { "dhA^- =Ak B T\a — = Avk B ThiN. (3) 
. L ha 
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FIG. 2. Examples of three dimensional figures without a 
length scale; grey surfaces indicate truncation only for graph- 
ical representations, (a) circular cone, (b) two-dimensional 
sector of a circle (in 3D space), (c) pyramid, (d) wedge, (e) 
plane, (f) line, (g) point. A polymer is to be attached to the 
point (g), to the apex of (a), (b), (c); or to any point on the 
edge of (d), or the entirety of (e) or (f). The plane and line 
can also be semi-infinite with a polymer attached to their edge 
or end-point. 



This work can also be computed from the change in free 
energies between the final and initial states, due to the 
change in entropy, as 



-TAS = TS C - TS C 



fcsT(7 c - 7cp )lniV, (4) 



with the entropy S = — fee In A/" computed from Eq. ([2]). 
By equating W and AJ- we find 



A 7c - 7cp 

A = = rj cp - 7} c ; 



(5) 



the final result obtained from the scaling law 

7a = (2 - V S >, (6) 

where the exponent r\ characterizes the anomalous decay 
of correlations (~ l/r d ~ 2+n ). 

The arguments presented in the previous paragraph 
rely only on the fact that configurations of the obstacles 
lack a length scale for both h — > and h — >■ oo. Conse- 
quently, similar reasoning can be applied to a variety of 
objects such as those depicted in Fig. [2] or combinations 
of such objects. The force prefactor is then related to the 
exponents in the initial and the final states by 



A = 



7 



initial 



7 



Una! 



= v 



„ final 



initial 



n 



(7) 



Consider, for instance, an ideal polymer in free space (or 
held by a point-like object (Fig.[2g)). Its number of con- 
figurations is z N , i.e. 7 = 70 = 1, and correspondingly 
r/o = 0, since v = vq = 1/2. Subscript will henceforth 
indicate exponents describing ideal polymers. If the end 
of the polymer is brought into contact with a plane, then 
in any space dimension d the number of configurations 
scales as z N N~ 1 ^ 2 , i.e. with 7 = 7 p0 = 1/2, and corre- 
sponding to ?7po = 1. Thus for a long polymer held by 
a point-like object at a distance h from a plane, the en- 
tropic force in Eq. ([1} will have a prefactor A = 1. This 



result is valid for any d. For self-avoiding polymers in 
d = 2, the exponents are 7 = 43/32, 7p = 61/64 H3 
and i/ = 3/4 Jill leading to A = 25/48 « 0.52, while in 
d = 3, we get 7 1.158 [2^, 7p » 0.70 and v « 0.59 (H 
leading to A ~ 0.78. Similar exact results can be derived 
for more complex self-avoiding polymers, such as star 
polymers or branched polymers of arbitrary topology in 
d=2[Hl. 



Note that the force amplitude does depend on the sur- 
face to which the polymer is attached. For example, con- 
sider the configuration in Fig. [TJi where the polymer is 
attached to the plate. The initial configuration corre- 
sponds to a polymer attached to an infinite surface, thus 



'1 



initial 



?7p, while the final configuration is the same as 
in Fig. [TJ;, with jy final = jy cp . This leads to A — rj cp — r) p , 
as opposed to A = r/ cp — r/ c for the case that the polymer 
is attached to the cone. There are a number of detailed 
studies of small surfaces (such as AFM tips) (3lTl37| com- 
pressing a single polymer (or few polymers) grafted to a 
surface. However, all these studies concern non-scale in- 
variant geometries. 

When the obstacles are separated by h, the loss of 
polymer entropy (and the corresponding pressure lead- 
ing to the force in Eq. (Q])) is concentrated in the con- 
finement region. The part of the polymer that wanders 
away from this area is relatively unperturbed and does 
not contribute to the force. Therefore for large N, Eq. ((T|) 
is independent of polymer size. This is not the case when 
the entire polymer contributes to the force, as for a poly- 
mer held between two plates, where the force is propor- 
tional to N. Thus, the argument fails when dimension- 
ality of the system is changed between the initial and 
final states. For example when a 3-dimcnsional polymer 
connected to a two-dimensional plane approaches a par- 
allel plane, in the final configuration the confined poly- 
mer is essentially two-dimensional. Free energies in the 
initial and final state have different extensive parts, and 
the polymer-mediated force depends on the number of 
monomers. 



It is worth reiterating the idealizations leading to 
Eqs. H]) and ([7]). By focusing on entropy, we have as- 
sumed that the only interaction between the polymer 
and the obstacles is due to hard-core exclusion. Attrac- 
tive interactions between the polymer and surface will 
introduce temperature dependent corrections, and addi- 
tional size scales. Weak interactions are asymptotically 
irrelevant, but strong interactions may lead to a phase in 
which the polymer is absorbed onto the obstacles, ren- 
dering the entropic considerations presented here inap- 
propriate. Real obstacles are self-similar over a range of 
length scales; for example an AFM tip may be rounded, 
or abruptly change its angle. Equation |T]) will then be 
applicable to a corresponding interval, and there ma y b e 
non-trivial crossovers between the different regimes [38| . 
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III. IDEAL POLYMERS NEAR 
SCALE-INVARIANT SURFACES 

Having reduced the computation of the force to calcu- 
lation of correlation functions, the remainder of the paper 
is dedicated to computing exponents characterizing poly- 
mer correlations. The behavior of idealized polymers, i.e. 
in the absence of self-avoidance and other interactions, is 
closely related to the diffusion problem: the spatial con- 
figuration of a polymer of iV-segments can be regarded 
as the space-time trajectory of a random walker viewed 
at integer times up to t = N. Summing over all trajec- 
tories starting at r and ending at r' is then equivalent 
to considering the probability of diffusion from r to r' in 
time t, which satisfies the diffusion equation 



dP(r,r',t) a , , N 

g t ' =DV 2 P(r,r',t). 



(8) 



The prime sign on the Laplacian indicates the spatial 
derivatives with respect to components of r', and the ini- 
tial condition is P(r, r', t = 0) = S d (r — r'). The diffusion 
constant D is chosen such that in free space the mean 
squared distance coincides with the random walk value 
of ((r — r') 2 ) = a 2 N = 2dDt in d space dimensions, and 
thus D = a 2 / 2d. If a polymer is confined by repulsive 
walls, its configurations correspond to paths of a random 
walker that does no cross these boundaries. In the lan- 
guage of diffusion, the exclusion of paths that cross the 
boundaries is accomplished by assuming that they are 
absorbing surfaces. 

To compute variations in the number of states, we 
are interested in the total survival probability S(r, t) = 
J P(r, r', t)d d r' . which also satisfies the equation j39j 



dS(r,t) 
dt 



DV 2 S(r,t), 



(9) 



with absorbing boundary conditions. The initial condi- 
tion for survival probability is S(r,t = 0) = 1, every- 
where inside the space where the particle can diffuse, 
and S(r,t) — on the absorbing boundaries. (While 
the space and time variables are continuous in the above 
equation, appropriately discretized equations can be ap- 
plied to diffusion or random walk on a lattice.) 

In the absence of confinement, the number of configu- 
rations J\f of an ideal polymer grows exponentially with 
its length N. On a regular discrete lattice, N is sim- 
ply related to lattice coordination number z by M = z N . 
The presence of repulsive boundaries reduces the number 
of configurations to 



Af(r,N) = z N S(r,kN) 



(10) 



where the prefactor k = a 2 /(2dD) converts the dimen- 
sionless number of steps N to time t in the diffusion 
equation. Numerical solutions of Eqs. (|8I9[) can be easily 
found for many geometries, and analytical solutions are 
available for a number of simple shapes [40]. The gen- 
erally complicated solutions simplify to some extent for 



the self-similar shapes depicted in Fig. [2] In a coordinate 
system centered on a special point of the object (such as 
tip of a cone) the dimensionless function S can only de- 
pend on the dimensionless vector w = vj\fDt. Thus, 
S(r,t) = H(w), and Eq. ^ reduces to 



V^tf + - w V W H = 0, 



(11) 



where the subscript w indicates derivatives with respect 
to components of w. In terms of these dimensionless 
variables, the function H vanishes on the absorbing sur- 
faces and increases to unity away from the surface (at 
separations of order one). For some geometries, the so- 
lution to Eq. (jlip can be presented in terms of a ra- 
dial distance w, and a combination of angular variables. 
For example, generalized cones in d-dimensional space 
are described by a single polar angle 6 and d — 2 az- 
imuthal angles <f>,ip,' ■■ . (While d-dimensional wedges 
can be treated similarly, it is more convenient to note 
that the solution is independent of the coordinates par- 
allel to the wedge, while the cross-section perpendicu- 
lar to the wedge resembles a cone.) Furthermore, for 
wCl the distance dependence is expected to be a sim- 
ple power law w Tl0} i>(d, (/>,...). In this limit, the second 
term in Eq. (|lip becomes negligible, and the problem re- 
duces to solving the Laplace equation V 2 v (w Vo< S) = 0. 
The boundary conditions are now implemented solely in 
terms of angular variables, and the problem is reduced 
to finding a non-negative function 'J that satisfies these 
conditions, and the corresponding i]q for which such so- 
lution exists. The dependence of the latter value on the 
boundary shape will be emphasized by denoting it as 7j s q. 
Since, S ~ w ns " = (r / \fl5i) ria0 for a random walk that 
starts a short distance |r| = a away from the cone, the 
survival probability is S ~ (a/ \f~Di) Va0 , and consequently 
the number of configurations of an ideal polymer grows 
as z N N-^°/ 2 = z N N^°- 1 . We see that 7s0 = 1 - i hQ /2 
as in Eq. © for vq — 1/2. 



IV. IDEAL POLYMERS NEAR CONES AND 
PLATES 

As an example, consider a polymer attached to the 
apex of a cone, as in Fig. QJ,, or to a point of contact 
between a cone and a plane, as in Fig.[TJ;. When general- 
ized to d-dimensions, this system has cylindrical symme- 
try, and ^ is independent of the d — 2 azimuthal angles. 
Therefore, the Laplace equation simplifies to 



1 



(sin6») d - 2 dO 



(sin 6>) 



d-2 



d^ 

~d£ 



+ VQ (d-2 + V0 )^(e) = o, 

(12) 



with an appropriate boundary condition on ^(0). For an 
isolated cone, the function ^> must be positive and regu- 
lar outside the cone, with d^/d6\g =v = to avoid a cusp 
on the symmetry axis, and 'I'(O) = on the cone surface. 
For the cone+plate, the appropriate solution is positive 
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FIG. 3. (Color online) The exponent r; s o for ideal polymers 
in d —2 (dot-dashed), 3 (solid), 4 (dashed) for cone ("s=c") 
of angle 6 (bottom curves), and "s=cp" (top curves). (From 
Ref. @.) 



and vanishes both at 9 — 9 and 9 — ir/2. The first case 
was considered by Ben-Naim and Krapivsky [4l| in con- 
nection with diffusion near an absorbing boundary [42j . 
and we follow their derivations. The solution in general 
d requires the use of associated Legendre functions, but 
simplifies in a few cases described below. 
• For d = 2, the problem of a cone coincides with that of 
a wedge. In higher space dimensions the wedge remains 
equivalent to the two-dimensional the function 

H is independent of the d — 2 coordinates parallel to the 
wedge. Thus, the following results for d = 2 are also ap- 
plicable to wedges in any d. Equation (|12p now reduces 
to ty" + r]o 2l & = 0; which is solved by linear combina- 
tions of sin(?7o6l) and cos(?7o0). The requirement that \P 
is positive and regular, and vanishes on the object(s), 
yields 



2(71-6)' 



and 



2tt 



26 



(13) 



Both results (depicted in Fig. [3j go to a finite value as 
O — > 0, reflecting the strong reduction in configurations 
due to the remnant (barrier) line, and 770 — > 00 when the 
boundaries confine the polymer to a vanishing sector. 

This computation can be easily generalized to the case 
of two cones (wedges) ("s=cc") with apex semi-angles 9i 
and 9 2 touching at their tips with a common symmetry 
axis. If one of the wedges is tilted by an angle t relative 
to that axis, then the survival probability S will be de- 
termined by the larger sector of free space. The resulting 
exponent is 



VccO 



7r-(ei + e 2 -r) 



(14) 



By proper choices of Gs and r we can reproduce the 
results of Eq. ([TB")) as particular cases. 



• For d = 4, the substitution ^ = u/s'm9 simplifies 
Eq. (|12p to u" + (?7o + l) 2 u = 0, solved by a linear com- 
bination of sin[(?7o + 1)#] and cos[(?7o + 1)0], and we find 



9 



9' 



and 7], 



29 



cpO 



29 



(15) 



depicted by the bottom and top dashed lines in Fig. [3] 
The cone exponent r] c0 vanishes linearly with 9 - a needle 
in four dimensions is "invisible" to a random walker. 

As in the two-dimensional case, we may generalize to 
two touching cones aligned along common axis with apex 
angles 9i and Q2, with 



»7cc0 



9i + 9 2 

-(9i + e a )' 



(16) 



Unfortunately, when one of the cones is tilted, the prob- 
lem loses its azimuthal symmetry and Eq. (|12p with its 
simple solutions is no longer valid. 
• For d = 3, Eq. (IT2l becomes 



d 2 ^ 



(1 - n z )— -2n—+ 770(770 + 1)* = 0, 



dfi 2 



(17) 



where fi = cos 9. The general solution to this equation is 
given by regular (rather than associated) Legendre func- 
tions 



V(e) = aiP, l0 (p)+oaQr l0 (jj). 



(18) 



Note that P a (l) = 1 f° r an Y a , while P a (—1) diverges 
for noninteger a. Similarly, Q a (±l) is divergent. The 
linear combination in Eq. (|18[) can be made regular at -1 
by a proper choice of a^fai- For the geometry described 
in Fig. QJi the solution must be regular for 9 < < it. 
Instead of using combinations of P and Q, we can simply 
use Pj7o(— cos 9), that will be regular at cos# = — 1. The 
value of 770 is then determined by requiring 



P„ (-cos9) = 0. 



(19) 



Since ^ cannot change sign in the physically permit- 
ted region, the smallest possible 770 must be chosen 
The corresponding solution has been described in de 
tail in Ref. 



41] 



In particular, in the limit 9 — > 0, 
the exponent 770 vanishes through a logarithmic singu- 
larity (~ 1/1 In @ I ) , while for 9 —> it it diverges as 
77 w 2.405/(7r- 0). 

For the geometry of a polymer attached to the apex 
of a cone touching an infinite plane (Fig. [1}:), we need 
to solve Eq. (fTT|) with *5>(9) vanishing at both 9 and at 
7t/2. The latter condition corresponds to /1 = and can 
be assured by setting in Eq. ([T8"|) the ratio 



(20) 



02 = P vo {0) = 2 

a i Qvo(Q) 7rtan(7r?7o/2) 



With the above choice, the angular function can be con- 
veniently written as 



#(0) = 7rsin(7T77 O /2)P ?;o (cos0) + 2 cos(7T77 /2)Q^ 



(21) 
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For non-integer a, P a (p) and P n (—ij) are linearly in- 
dependent and both solve Eq. p7|) [43j ]. Thus for non- 
integer 770, Eq. (|2"Tj) can be replaced by 



*(9) = P Vo (caa9)-P Vo (-caB 



(22) 



Again, we must choose the smallest 770 f° r which this 
function vanishes for 9 = 0. The resulting exponents 
(which cannot be cast as simple functions), are plotted 
as solid lines in Fig. |3l From bounds on the roots of 
Legendre polynomials piHia ]. we find that in the limit 
— > 7r/2, the exponent diverges as 770 ~ 2tt/(tt — 26), 
while as — > 0, 770 — > 1- Indeed, the solution in Eq. (f2"Tj) 
has no roots for 770 < 1. Asymptotic expressions for 
Pa(fJ') and Q Q (/x) near \x = 1 are known and can be used 
to determine positions of roots in such limit, indicating 
that 770 w 1 + 1/1 ln(G)|. 

It is worth noting that in all cases above (d = 2,3,4) 
the cone+plate exponent has the identical divergence 
770 ~ 2tt/ (tt — 20) as — > ir/2. This may be justified by 
arguing that a diffuser confined between an almost flat 
cone and the plate encounters absorbing boundaries not 
dissimilar to that of a wedge. 

• In general d, with the substitution ^£(9) = sin -5 9 u(/z), 
where /1 = cos 9 and d = (d — 3)/2, Eq. (TT2"j) becomes 



dpJ 



2/j 



du 
dfi 



(t]o + 5)(tio + 5 + I) - 



1 



u = 0, 



(23) 



which is solved by associated Legendre functions [43| 
that can be chosen in several forms such as P^ s +S (±fj,), 

0^+ or their combinations. The particular choice 
of one or pair of linearly independent functions depends 
on the boundary conditions and the specific values of 770 
and S. For the case of cone in different dimensions the 
reader is referred to Ref. [4l|- To find a proper form for 
u in (integer or fractional) d dimensions, we can closely 
follow the case of d = 3: For a single cone we may chose 



P 



-8 

Vo+S y 



-cos#), since sin - 9 P~ o+s (— cos9) has no 
cusp at 9 = tt. The value of 770 will be set by the re- 
quirement that the function vanishes for 9 = 0. For 
cone+plate configurations we may choose 



(24) 



which for a\ = a 2 vanishes on the plate {9 = tt), and will 
be a suitable solution as long as 770 is not integer and is 
chosen such that the function vanishes for 9 = 0. 

For all d, the cone becomes a plate for = tt/2. 
Correlations with one point approaching a surface are 
easily obtained by the method of images [It} leading to 
Vco = Vpo = I7 which is clearly seen from the intersection 
of the curves in Fig. [3] 

Note that in d = 3 both cone and cone+plate expo- 
nents approach their limit — > via a logarithmic sin- 
gularity, while for d — 4 they approach that limit linearly. 



For intermediate dimensions 3 < d < 4, the limiting be- 
havior as — > is given by 



VcO 



VcpO 



r(l-e/2) 



1 - 



1 



^rr(l/2-e/2) 

4T(2-e/2) 
+ v^r(l/2-e/2) 



and 



1 - 



(25) 



where d = 4 — e. The singular behavior in the above 
equations is a power law pt> with pa = d — 3, which 
has a simple geometric interpretation. The ideal poly- 
mer can be regarded as a self-similar object with fractal 
dimension 2, while the remnant of the cone as — > is 
a (semi-infinite) 1-dimensional line. When embedded in 
d-dimensional space, the intersection of the two entities 
(random walk and remnant line) spans a space of dimen- 
sion d— (2 + 1) = pq. Indeed, for d < 3 the limiting value 
is different from the case without any cone indicating the 
finite probability of intersection of the polymer with the 
semi-infinite barrier line. 

In both d = 2 and d = 4, we found simple expres- 
sions for the ideal polymer attached to the contact point 
between two coaxial cones of opening angles 0i and 
02. However, in general d the solution has the form in 
Eq. (J2U). The vanishing of u on the first cone requires 
ai/a,2 = cos(0i)/ cos(— 0i). With this choice, the only 
free parameter left is the exponent 770. To satisfy the 
condition u{9 = ©2) = 0, we must have 



P Vo+S 1 



cos0 1 )P-^(cos0 2 )=O. (26) 



Unfortunately, values of 770 satisfying the above equation 
cannot be expressed as simple functions. However, for 
small 0i and ©2, we can prove that 



7?cco(©i, ©2) ~ ?7co(@i) + ?7co(©2), 



(27) 



i.e. the exponent for two sharp cones is approximately 
the sum of exponents for individual cones of opening an- 
gles @i and 2 . 



V. SELF-AVOIDING POLYMERS: 
SIMULATIONS 

Self-avoiding walks (SAWs) provide a convenient 
model for exploring universal aspects of swollen (coil) 
polymers with short-range interactions. The number of 
SAWs in d = 3 dimensions space (without obstacles) is 
governed by the exponent 7 w 1.158 [29( (77 0.03), 
while in d = 2, 7 = 43/32 48] (77 = 5/24). (For ideal 
polymers 70 = 1 (770 = 0) at any d.) There are a number 
of results regarding 7 S (with obstacles) for polymers con- 
fined by wedges in 2D and 3D [13, Hfl \M^M- A SAW 
anchored at the origin and confined to a solid wedge (in 
3D) or a planar wedge (in 2D) has an angle-dependent 
exponent 7wedge that diverges as the confining angle van- 
ishes. Studies of a polymer attached to the tip of a 2D 
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sector in 3D, and to the apex of a cone were reported 
in Ref. [53l . Extensive analytical [H], and numeri- 
cal [m Ell studies of SAWs anchored to a solid plate in 
3D find 7 S = 7 P in the range of 0.70 [IH to 0.68 |i] 
(r/ p = 0.81 or 0.84, respectively). (For ideal polymers 
7p0 = 1/2 (?7po = 1) for any d.) We are not aware of 
specific results pertaining to the geometries depicted in 
Fig. [TJ and provide below our numerical and analytical 
estimates. 

We performed numerical simulations for SAWs of 
lengths N = 16, 32, . . . , 1024 on a cubic lattice. The 
SAWs were generated by a dimerization method [U HH 
in which an unbiased TV-step SAW is created by at- 
tempting to join two TV/2-step SAWs previously obtained 
by the same algorithm: If the resulting A-step walk 
does not intersect itself, the new configuration is ac- 
cepted. If the combined walk has intersections it is re- 
jected; both TV/2-step walks are discarded and the pro- 
cess is repeated. As the rejection rate increases slowly 
with TV, this is an extremely efficient method for cre- 
ating unbiased collection of SAWs. We generated 10 8 
SAWs for different values of N; each SAW is attached 
to the origin and checked whether it touches a cone (or 
cone+plate) . The probability that an TV-step SAW does 
not intersect the excluded space is the ratio of the per- 
mitted number of walks to the total number of SAWs, 
i.e. p N = M s /M ~ N^-^/N^- 1 = N A ^. The ratio 
Pn/p2N — 2 A7s is then used to estimate A7 S for each 
sequential pair of As. The results are extrapolated to 
TV — » oo by plotting the estimates versus 1 / y/N ; errors 
are due to both finite sample size and insufficiently large 
A. The estimated exponents A?7 S = r/ s — rj = (7 — 7 s )/f 
(full symbols in Fig. 2]) are somewhat lower than the cor- 
responding values for ideal polymers shown by dotted 
lines. The difference primarily originates in the change 
in r] s as a result of self-avoiding interactions. (It is slightly 
moderated by a small shift of rj in free space.) Figure 3 
of Ref. [23| presents A7 S as a function of 0, in which 
form results for SAWs are much closer to those of ideal 
polymers due to multiplication by the self-avoiding and 
ideal values of v, respectively. 

The dashed line in Fig. @] depicts the result [53| of an 
e = 4 — d expansion in which the 2-dimensional surface 
of the cone is treated as a weakly repulsive potential. 
The lowest order result, Arj c = (3e/4)sinO, captures 
some features of the numerics for < tt/2, but fails 
dismally for > 7r/2 (since the polymer simply jumps 
to the larger space inside the cone). It also fails to cap- 
ture the correct behavior as — > 0. A better approach 
is to exclude the interior of the cone, as initiated by 
Cardy [13, H(| for the wedge geometry. The correspond- 
ing computations for a cone, described next, are more 
complicated, and build upon more recent results pertain- 
ing to the electrodynamic Casimir force on a cone [59j . 



VI. SELF-AVOIDING POLYMERS: 




FIG. 4. (Color online) Dependence of the exponent differ- 
ence Arjs = r/s — n on angle O for the cone+plate (top curves) , 
and an isolated cone (bottom curves). The dotted curves are 
the exact values for an ideal polymer, while the data points 
represent numerical results for SAWs in the same geometry. 
All error bars are estimates of the uncertainty in extrapola- 
tion to TV — > 00. The dashed line depicts prediction of an 
e-expansion in Ref. [k| (with e = 1) for 'weak surface repul- 
sion.' The solid lines show the values of r) c from Eq. (|38|) . and 
?7 C p from Eq. (|50|) . with e — 1. For the latter, Arjs = rj s since 
n = + 0(e 2 ) in free space. 



EPSILON EXPANSION 
A. Cone 

The constraint of self-avoidance is relevant in dimen- 
sions d < 4, where exponents are computed perturba- 
tively in e = 4 — d. A primary element of such calculation 
is the Green's function in the absence of interactions (self- 
avoidance), but in the presence of obstacles. For this we 
need the full solution to the Laplace equation in four di- 
mensions, going beyond the large separation asymptotic 
form obtained via Eq. (|12p . (For a general e-expansion, 
one has to know the Green's function in 4 — e dimensions. 
However, we carry out this expansion only to first order 
in e. Since the strength of the interaction is proportional 
to e, it suffices to obtain the Green's function in d = 4.) 
Subject to appropriate boundary conditions, the Green's 
function is the solution to 

-V\ 4) G {x,x') = 8\x-x'), (28) 

where x and x' are spatial coordinates, and the subscript 
4 emphasizes the space dimension. Since we are inter- 
ested in conical boundaries, it is convenient to break up 
the Green's function along the polar angular coordinate. 
In accomplish this, we need a complete set of functions 
of the other three coordinates including the radius. One 
can exploit the Kontrovich-Lebedev transform [6(| to for- 
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mulate the completeness relation. The procedure is anal- 
ogous to Ref. [59| where some of the authors have carried 



out a similar analysis. Details of the calculation are given 
in the Appendix, and the final result for the Green's func- 
tion for a single cone is 



Go{x,x') =J2 



(-l)<7rr(p n + Z + l) 



„P"-! P 



-1-1/2 



(- cos 8)Y lv M, 0) P- l I${- cos 6')Y* m {f, 0') 



2sin(/9„7r)r(p n - I) g p~ 



-1/2 
-1/2 



(-cos 6) r 



+i 



Vsii 



Vsin^ 7 



(29) 



where the sum is over the triplet of integers n > 0, I > 0, 
and — / < m < +/. The important exponent p n is the 
n-th root of the transcendental equation 



P 



-1-1/2 



Pn-1/2 



(-cose) = o. 



(30) 



In the above equation, r < and r > refer to the smaller and 
larger radial coordinate for x and x' . Such division of the 
Green's function is appropriate since one endpoint of the 
polymer is close to the tip of the cone and the other end 
is far away. 

The interacting Green's function is obtained by sub- 
tracting polymer configurations that self-intersect. In the 
perturbative analysis, the strength of self-repulsion is in- 
dicated by a parameter u, and to lowest order one must 
subtract contributions forming a single intermediate loop 



I 

symmetry and yields 1(1 + 1) when summed over all spher- 
ical harmonic functions of degree I. In performing the 
integration over the radial coordinate r" , we need to ex- 
tract a logarithmic contribution which is ultimately ex- 
ponentiated to a power-law. Similar to the analysis in 
wedge geometry by Cardy, such a logarithm comes only 
from the region r < r" < r' . There remains an inte- 
gral over the polar angle 9" which is cumbersome since it 
requires finding the roots of the transcendental Eq. (|30[) . 

We shall not dwell further on the e-expansion for a 
general opening angle. Instead, we focus on the limit of 
a sharp cone as 6 — > 0, where the leading singularity 
comes only from the I — channel of the loop Green's 
function. This is because a sharp cone couples only to 
the lowest spherical partial wave. The loop propagator 
is then a sum over n: 



Gi = Go- u d*x"G (x iX ")G r (x",x")G (x",x'). (31) G r \ l=m=0 (x" , x") = 



1 sin 2 n7r^4- - sin 2 n9" 



4 7r 2 r //2 Z-^i n 



Sill 



In the language of quantum field theory, the first term is 
the "free" Green's function, albeit in the presence of ex- 
ternal boundary conditions, while the second term is the 
one-loop correction. Note that the intermediate Green's 
function is regularized by subtracting the Green's func- 
tion in empty space, Go, hence the superscript r, i.e. 
G = Go — Go- In a renormalization group treatment, the <x 
parameter u changes with scale of observation, ultim ately 
arriving at a fixed point value of u* — 2n 2 e [27], H3, l53j . 
This universal value is a characteristic of the polymer and 
independent of obstacles or boundary conditions. To cal- 
culate the scaling properties for r 3> a, it is sufficient to 
include only the first term (n = 1, I = m = 0) from 
Eq. (|2"5]) in the non-loop propagators 



(33) 



Equations (|3"2"|) and (f3"3"|) are used to perform the integral 
over 9" and summation over all n. The radial dependence 
of the Green's function is then obtained as 



»?rf>+2 



1 + eln— f — 61n6 

r> \4-7r 



.166 



, (34) 



where the term proportional to e in the bracket is the 
one-loop correction to the Green's function. In d = 4 — e 
dimensions, t] c q is given by Eq. (|25[) . Next we exponen- 
tiate the radial logarithm to obtain the renormalized ex- 
ponent, 



Go\i,m=o( x > x ') 



sinpi (tt — 6) sinpi (tt — 9') 



Vc = VcO + 



91n9 

47T 



.166 



(35) 



47T 2 ' 



-P1+1 



(32) 



(As in the wedge computation by Cardy [271 . |50| , higher- 
order terms give corrections in higher powers of a/r.) 
The intermediate loop, however, can be of any size re- 
quiring the entire sum (albeit regularized by subtracting 
the result for 6 = to remove an unrelated divergence). 

The loop correction is an integral over the whole space: 
The integration over angles <t>" and ip" is trivial due to 



As the first correction to rj in empty space appears at 
the order of e 2 , the above (first-order) result for the cone 
vanishes logarithmically as 6 —> 0. As in the case of the 
radial logarithms in Eq. (|34[) . it is tempting to interpret 
the angular logarithm in Eq. (|35l) as a signature of a 
power-law for 6 — > 0. Expanding r) c0 to the first order in 



6 

VcO ~ — 

7T 



6 In 6 



- .22 6 e 



(36) 



9 



we find another In Q, originating from the expansion in 
4 — e dimensions of the phantom polymer (as opposed to 
the perturbative terms from the one- loop computation). 
Summing both contributions, we obtain 

Vc= ~ U- jelne-.06e) , (37) 
which can be recast as a power law 

?7 C (~ - .06 e 1 "^ . (38) 

As expressed in Eq. (f^5|) . the exponent ?7 c o vanishes with 
cone angle as O po with po = d — 3 = 1 — e for the phantom 
polymer (Fig. 13]). The above equation indicates that the 
vanishing of rj c for a self- avoiding polymer is governed 
by the modified exponent p = 1 — 3e/4. This depen- 
dence is shown for e = 1 as the lower solid line in Fig. [4] 
We may interpret this result as follows: self-avoidance 



swells the polymer, reducing its fractal dimension from 2 
to^~ 1= 2 — e/4 to lowest order. The dimensionality of 
the intersection of the fractal coil with the remnant line 
is now d — 1 — v~ =1 — 36/4 = p. It would be inter- 
esting to see if this connection holds in higher orders in 
perturbation theory. The current numerical results are 
not accurate enough to test this conjecture. 

Finally, we note that for two touching cones of small 
angles Oi and 82, we find the analog of Eq. (|2"T1) to first 
order in the e-expansion, namely 

??cc(e 1 ,e 2 )«r ?c (e 1 ) + r ?c (e 2 ). (39) 



B. Cone-Plate 

When constrained by a cone touching a plate, the 
Green's function is given by 



(^-ffi-cosfl) - P;*Zl',l(cm6))Y lm ty,<l>) (p-'lffl-cosfl') - P^rfflrasg'))^', <j>>) 
V sin 9 V sin 9' 



(40) 



where the sum is again over the triplet of integers n > 0, I > 0, and —l<m< +1. The exponent p n is now the rt-th 
root of the transcendental equation 

(" cos ) - Ppl-l/t ( cos ®) = 0- (41) 
Again for the non-loop propagators, it suffices to keep only the first term in the sum 

n l»=i r- ->\ - r <~ 1 sinpiWS-^sinp!^-^) 



> 



However, the intermediate loop should be summed entirely. Here again, we focus on a sharp cone. It is more convenient 
to express the regularized Green's function for the intermediate loop as 

G 7 = Gq — Gq = (Gq — Gq) + (Gq — Go) , (43) 

where Gq is the Green's function in the presence of the plate alone and is independent of the cone angle 6. Also note 
that Go is the Green's function in empty space. Each bracket on the RHS of Eq. (|43p can be treated separately. The 
first bracket, to the leading order for a sharp cone, is given by 

oo • 2 ir/2 — 8" • 2n all 

1 ^2, 1 sin mr ,„ a sin 2n9 

(Go - G%) \ M (x",x") = ' J2 Z ^T^, ■ ( 44 ) 

n—l 

Also the second bracket in Eq. (|4"3")l is merely the Green's function in the presence of a plate. The latter can be cast 
as a series expansion by using the method of images: 

(Gg - Go) = -E ^n+i ^rfa P "-^ /2( ~ ^ S e)p-!_-ji 2 (cos9') Y lm {1>, m%W, <P')- (45) 

nlm > 



I 

The Green's functions appropriate to the external legs the limit of O — > 0, 
and the loop should be inserted in Eq. (|3Tj) to obtain, in r; cp0 

Gi (X -77 

r l)cp0+2 



l + eln^r-i + ^ + .66e 



r 



> 



(46) 
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where the exponent 77 cp o is given by Eq. (|25[) . Upon ex- 
ponentiation in r, we find 



Vcp = VcpO + 



i eine 



.66 6 



(47) 



Note that the loop-correction does not vanish for 6^0, 
instead going to the limiting value of — e/8 due to the 
presence of the plate. Expanding 77 cp o to the first order 
in e, 

7T V I" / 

we find a logarithmic dependence in O. Summing both 
contributions, the exponent r\ to the lowest order in e and 
becomes 



46 



?7c P = 1- - + — 1 - T eln0- 

O 7T 



3 



(49) 



As in the previous section, we exponentiate the logarithm 
in G to obtain the power-law dependence 

77 C p«l-^+ (^-Me^je 1 -^. (50) 

Revealingly, the exponent p = 1 — 3e/4 is the same as 
in the case of a single cone, whereas the amplitude has 
changed. Furthermore, it approaches a constant value 
in the limit of vanishing angle. This result is shown for 
e = 1 as the upper solid line in Fig. 2] 

Finally we can compute the amplitude of the force ac- 
cording to Eq. (O, using Eqs. (1581 and ((501 . as 



A = l- - + [ -- .80 e ) Q 1 -*', 

8 V 7T 



(51) 



C. Star Polymers 

From a practical point of view, the amplitude A = 
Vcp ~ Vc is typically a number of order unity. Thus the 
force in Eq. (fTJ) at room temperature is roughly around 
0.1 pN at 0.1 fim separation; at the margins of possi- 
ble measurements by current force apparatus. The force 
can be increased by attaching several polymers to the 
tip. For ideal (phantom) polymers the force is enhanced 
by /, the number of polymers, while interactions modify 
this conclusion. The interactions are of two kinds: self- 
interaction of a single arm (intra-arm) and interactions 
between two different arms (inter-arm). The former con- 
tribution is the same as computed above to first order 
in e, and should only be multiplied by / for an /-arm 
polymer. The latter, however, gives rise to a new term 
multiplied by /(/ — l)/2, the number of pairs. Thus to 



first order me, r/ — fvo+fv 



i i f(f-i)„e 



rf , where r/o denotes 



the exponent of a single phantom polymer, and r/ 1 and rf 
are contributions from intra-arm and inter-arm interac- 
tions respectively. In the absence of a cone, this set-up is 
similar to the widely studied case of star polymers [61] . 



To zeroth order (i.e. in the absence of interactions) the 
correlation function describing a pair of polymers start- 
ing at the same point x, but with different endpoints at 
x'i and x' 2 , is simply the product of two free Green's func- 
tions, Go <§) Go ■ Subtracting configurations in which the 
two polymers avoid each other yields to first order the 
interacting Green's function 

Gf Y = G G -u I dV'Go>,a/0Go(Azi)Go(zV2), 



(52) 

where we have summed over the intersection point x" . 
The prime on the Green's function indicates that we 
have only considered the inter-arm interactions. Also 
note that no regularization is needed since the integral 
in Eq. (|52|) is finite. For a long polymer, it is sufficient 
to include only the first term in the Green's function, i.e. 
Eq. l|3"2l for a single cone and Eq. (|4"2l for a cone-plate 
configuration. The leading logarithmic contribution for 
a cone is then 



G(2) >o 
i — <jo 



1 + em- 



27 



(53) 



where r max depends only on r[ and r' 2 . Hence, by expo- 
nentiating the radial function, rf for a cone is obtained 
as 



Vc= e 



1 



2^ 



(54) 



The corresponding e-expansion for cone+plate gives 

79 



Vcp 



3tt 



(55) 



These equations then lead to the force amplitude 



AS) 
f 



, 80 + il(f-l)U e 1 - 36 / 4 



12tt 



(/-I) 



(56) 

Interestingly, interactions amongst the polymers reduce 
the force amplitude per polymer. This equation is similar 
to Eq. (I51|) with the addition of inter-arm interactions. 

Note that the same exponent p = 1 — 3e/4 again gov- 
erns the approach to a finite limit in Eq. (|55|) as & — » 0. 
Yet another experimental set-up is that of a polymer at- 
tached to the tip of one cone approaching a second cone. 
The entropic force is of course much smaller in this case, 
and we have verified that it also vanishes with cone angle 
as O p , bolstering the conjecture that this is a universal 
exponent related to the needle geometry. 



VII. DISCUSSION 

In summary, we have demonstrated that polymers ex- 
ert an entropic force AksT/h on a cone tip, with a 'uni- 
versal' amplitude A dependent on geometry, interactions, 
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and internal topology of the polymer. We conjecture that 
the singular behavior of the amplitude on vanishing cone 
angle is described by a new exponent, simply related to 
the fractal dimension of the polymer. There are many 
other self-similar shapes where a similar force law is ex- 
pected on the basis of scaling at length scales shorter 
than an appropriate correlation length. 

In this work we concentrated on situations where the 
polymer is attached to a single surface such as described 
by Fig. [lb or [TJl. However, our results can be easily ap- 
plied to a variety of other situations. First, we note that 
our approach is easily generalized to a slightly more com- 
plicated situation with one polymer attached to a cone 
while another is attached to a surface: The resulting force 
coefficient A is given by Eq. ([5]) with 7y lmtlal consisting of 
the sum of r] c for the two cones, while ?y final corresponds 
to two polymers (or a star polymer with / = 2) attached 
to the contact point between a cone and a plate ( "s=cp" ). 
An even more relevant situation is the force constant of 
a single polymer attached at both ends to a cone and a 
plate as in Fig. [T£. One may view this situation as an 
elaboration of the previous case of two polymers, when 
their free ends are joined to each other. Such a modi- 
fication significantly changes the behavior of the system 
when h is of order or larger than Rq. In fact the force 
changes sign at h ~ i?o, when the polymer state changes 
from stretched to compressed. For h <C Rq connecting 
the polymer ends modifies the free energy of the system. 
However, it is plausible that in such limit the force trans- 
mitted by the polymer (s) is not influenced by end-point 
connection, and therefore the force amplitude may be the 
same as the case of two disconnected polymers. Verifi- 
cation of this point requires a further study of the forces 



between the polymer and the boundaries. 
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Appendix A: Green's functions 

We first derive the Green's function for the Helmholtz 
equation in four dimensions; the corresponding function 
for the Laplace equation is then obtained as a limit. The 
former satisfies the Helmholtz equation with a delta func- 
tion source: 

- (V( 4) + fc 2 ) G (x, x') = S 4 (x - x 1 ) ; (Al) 

the latter is obtained by setting the wavevector k to zero. 
Consider spherical coordinates x = (r, 9, ip, 4>), where r is 
the radius, 9 is the polar angle and <f> and ip the remaining 
angular coordinates — a constant 9 specifies the surface of 
a cone in four dimensions. The cone is spherically sym- 
metric with respect to ip and <p (similar to the azimuthal 
symmetry of the three-dimensional cone). We expand 
the Green's function as 



G (x,x')=J2 At /■ !■ = 9im(r, e,r', 9') YU^m^', <!>')■ (A2) 



'rr' v»m sin f 

The prefactors are introduced for later convenience. With this ansatz, Eq. (IA1[) takes the form 

V (2) ~ {l t l ' 2 I + k ") 9im = fc^cosfl - cos0'). (A3) 
y ' r z sm9 z J r z 

Note that the Laplacian in the above equation acts in two dimensions, i.e. V? 2 n = ipid r r 2 d r + -^-gdg sin9dg. 

Given that the boundary is of conical shape, it is convenient to break up the Green's function in the coordinate 
9. We should then find a completeness relation in the function space of the other variable r. This is obtained by the 
Kontrovich-Lebedev transform 16011 



d\\smh\k iX - 1 /2(r)k lX _ 1 / 2 (r') = S(r - r'), 



where k v is the spherical Bessel function of order v. Using this relation, one can show that 

9im = I (-1)' J™ dXX fca-i/ 2 (*r)fc a -i/2 K) F ^ A ( + *_ + 1} i^lrfflcos B^P^i- cos fl>) , (A4) 

solves Eq. (|A3[) . Here, k is the imaginary frequency (k — in), 9 K = mm(9,9') and 6 > — max(9,6'). Note that 
this equation satisfies the Helmholtz equation in empty space. We discuss the boundary condition below. The 
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completeness relation, together with an identity regarding the Wronskian of Legendre functions, can be exploited to 
see that Eq. (IA4|) indeed solves Eq. (|A3|) . By a construction similar to Ref. (59|, we can analytically continue the 
complex order of the Bessel and Legendre functions to the real axis to obtain 



oo 



E, \, / nT(n + l + l) 1-1/2, a \D- l - 1 / 2 t a'\ 
«n-i/2(KT , <)fc„_ 1/2 (Kr > ) — _ — P n -i/2 (~ cos 6»)P_ n _/ /2 (- cos 6» ). (A5) 

n=l+l 1 [U L > 

Note that the asymmetry has shifted to the radial variable while the angular coordinates are treated symmetrically. 
We have used the symmetry P^-xi'i' (~ x ) = ( — ^-) l+n+1 P n -i/2 ( x ) ^ or integers I < n to restore the latter symmetry. 
Next we take the limit k — > to find 

' r(n + / 4- 1W"- 1 P^)' 2 (-cos8) P~ l 'U 2 (-cos9') 
G„(,,,') = EE E r(n n ^TT n - 1/2 r— ft - YU^)Y?mW,n (Afi) 

n=l i=0 ro=— I ? v v 

This equation provides the Green's function for the Laplace equation in empty space. For the Green's function in the 
presence of the conical obstacle, we modify Eq. (|A4J) to 

9im = J dXX fea-1/2 (Kr)fe iA -i/2 {Pi\-xjh™ s6 <) P i\-ij2(- cos0>)- 



p-'- 1 /, 2 (cose) 
^A-v^C-cose) 



p-^/^-cos^P^y^-cose')) , (A7) 



which vanishes when either angle is equal to G. It can be easily checked that this is indeed the Green's function: 
Upon acting with the Helmholtz operator, Eq. (| A3|) . the second line of the last equation gives a delta function while 
the third line vanishes since it has no discontinuity along 8 = 8'. Also the boundary condition is clearly satisfied 
as 8 — > O, the opening angle of the cone. Next we rotate to the real axis and take the limit n — > to obtain the 
corresponding Green's function for the Laplace equation 

,-1-1/2 



P " ', . = Y im W, <t>)Y? m W, n (A8) 

Vsint/' 

The integer n represents the n-th root of the transcendental equation 

p ;»- 1 1 / / 2 2 (-co S 9) = 0. (A9) 
The solution for p n also depends on and I. In the limit of 6 — > 0, i.e. for a sharp cone, it can be shown that [62[ 

Pn~l + n + - r ( 2 | + " + 1 ) ( tan 8)«+i (A10) 

F r(Z + 3/2)F(Z + l/2)F(n) V y V ' 

and the asymptotic behavior in this limit is dominated by the root for 1 = 0. To proceed further, we take advantage 
of the identities, 



P~ p % (cos 9) = -J — =- sin( P 0) , 
p 1 p V 7rsm 8 



the second of which dictates 



P;_Y; 2 (-co S e) = -1-^^-6)), (All) 



P„ = ^- (A12) 

7T — B 
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Hence, using Eqs. (|A8HA12|) . we find, for a single cone, 



Air z nrt sine/ sinfl' 

n > 

We can similarly find the Green's function in the space between two cones aligned along a common axis which touch 
at their tips. Designating the opening angles by 0i and 7r — 02, the available space is characterized by the polar 
angle 0! < 9 < 2 . The Green's function (for the Helmholtz equation) is similar to Eq. (|A4[) with the substitution 

of p &-i/2(™ <) p a-i/2(- cos6l >) h y 

D -i-X/2/ a \ P o!-i/2( cos@ i) -1-1/2, Q A {--1-1/2, „ , P ^!-i/2(- cos02 ) D -i-l/2/ fl n 

i/o (cos 0<) i—rk P\ i A (— cos0< P., ,;„(— cos6/>) , ^ /2 P.. ,/„ cos#> 

iA-X/2 V. <> P^Zlj 2 (~cos6i) »A-X/2 V <--/ J I zA-X/2 V >) P^j/^ (cos 6 2 ) »A-X/2 V >/ 



2 _ -iA-1/2 J %X-l/2 

p iA-i% 2 (- cos © i ) pr^/ficoB e 2 ) 



(A14) 

We then rotate to the real axis where we should find the roots of the transcendental equation 

P p- x"/2 /2 (- cos ©i)^ P - V2 /2 ( cos ©2) - "//'(cos 0i)P;_V// 2 (- cos 2 ) = 0. (A15) 

As we are mainly interested in a cone attached to a plate, we choose 0i = and 2 = ir/2. By rotating to the real 
axis and taking the limit k — > 0, the full Green's function is then obtained as 



Go(« *') = YY Y -r-D^^n+l+i) ^g(^e) 

° l ' ] ^ ^ ^ ^ ^ sin(p^)r(p„ - dpn (P+Ifr- cos©) - P-^CoobG)) 



n 1=0 m= — l 

p- l ~}£(-coa9) - P~ , ~ 1 / / „ 2 (cos0) P~ l ~^{- cos0') - P^^cos^) 

p„ — X/2 1 Pn — X/2 V ' p„— 1/2 V ' p ra — X/2 V ' ,, /, nvr* / ,/ ,/n / . -, „x 
" 7^ ; ' i Y lmWA) Y lmW A), ( A16 ) 

Vsmt/ vsmt/' 
where p„ is the n-th root of the transcendental equation 

p P T-i/l (- cos e ) - p P ; -x/2 2 ( cos e ) = °- ( Ai7 ) 

For a sharp cone, we can focus on I = 0, in which case, using Eq. (|A11|) . we find 

*-l£=e- < A18) 

The Green's function for the cone+plate configuration then reads 

r \ 1 \ ^ r <' 1 ' 1 ^Pn(^/2~e)sm Pn (TT/2~e') 

G \i =m=0 {x,x ) - — y — — — — . (A19) 

4tt z / 77, sin 6/ sin»' 
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